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Abstract 

We investigate in a model-independent way to what extent one can perturb tri-bimaximal 
mixing in order to generate a sizable value of |f/ e 3| ; while at the same time keeping 
solar neutrino mixing near its measured value, which is close to sin 2 #12 = 4. Three 
straightforward breaking mechanisms to generate \U e 3\ ~ 0.1 are considered. For charged 
lepton corrections, the suppression of a sizable contribution to sin 2 12 can be achieved 
if CP violation in neutrino oscillations is almost maximal. Generation of the indicated 
value of I t/ e 3 1 ~ 0.1 through renormalization group corrections requires the neutrinos to 
be quasi-degenerate in mass. The consistency with the allowed range of sin 2 #12 together 
with large running of \U e 3\ forces one of the Majorana phases to be close to it. This 
implies large cancellations in the effective Majorana mass governing neutrino-less double 
beta ((/3/3)oiy-)decay, constraining it to lie near its minimum allowed value of mo cos2#i2, 
where mo 0.1 eV. Finally, explicit breaking of the neutrino mass matrix in the inverted 
hierarchical and quasi-degenerate neutrino mass spectrum cases is similarly correlated 
with the (/3/3)oi/-decay effective Majorana mass, although to a lesser extent. The implied 
values for the atmospheric neutrino mixing angle 823 are given in all cases. 
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1 Introduction 



It is a remarkable achievement of experimental neutrino physics to have identified the leading 
form of lepton mixing, or Pontecorvo-Maki-Nakagawa-Sakata (PMNS), mixing matrix [1] U: 



U ~ Utbm P , where C/tbm 



V 



\ 
/ 



and P = diag(l, e lQ2 / 2 , e iCl3//2 ) contains the Majorana phases [2,3]. The above matrix Eq. (TjQ) 
defines tri-bimaximal mixing (TBM) [4]: 
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Currently, the values of sin 2 9 12 and sin 2 023 obtained from global fits of the neutrino oscillation 
data are indeed very close to those predicted by TBM [5]: 
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Here we have given the best-fit values as well as the la and 3a ranges (see also [6-8]). 
Obviously, even if Nature has chosen TBM as the lepton mixing scheme, one expects devi- 
ations from it on very general grounds. Straightforward examples are charged lepton correc- 
tions (i.e., corrections stemming from the diagonalization of the charged lepton mass matrix), 
renormalization effects, or explicit breaking in the neutrino mass matrix giving rise to TBM. 
Interestingly, in what regards the third mixing angle #13, a weak indication towards a non- 
vanishing value has recently emerged from a combination of two independent hints in solar, 
reactor and atmospheric neutrino data. Reference [5] quotes the following best-fit value and 
la range: 

" ? 13 = 0.016 ±0.010, (4) 



or \U e3 \ = sin 6»i3 = 0.126; 
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Vanishing #13 is thus disfavored at 1.6cx. 



Similar values and ranges have been found in other, independent analyses [10]. We note that 
the hint in the atmospheric data has been questioned [11], but that the recent MINOS data 
show an excess of electron events [12], which may be interpreted [13] as another hint for a 
non-zero 013. 

The allowed ranges for the mass-squared differences from the current global fit performed 
in [5] are 



Am| = 7.67_°;K:! x 10-5eV 2 



I Am\\ 



9 on+0.11, 0.42 10 _ 
z - oy -0.08,0.33 x iu 



(5) 



Note that the sign of |Am?| — (Ato 2 ^ ~ |Am 2 2 |, i.e., the ordering of neutrino masses is still 
not known. Regarding the neutrino mass scale, there are mainly three possibilities: normal 
hierarchy (NH) with m\ <C mi <C 7713, inverted hierarchy (IH) with 7713 <C m\ ~ mi, or 



1 The experimental results are so close to TBM that parameterizations of the PMNS matrix with TBM as 
the starting point have been proposed [9]. 
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quasi-degenerate neutrinos (QD) with m 2 , = m| ~ m 2 — m§ 3> Am|),|Am^|. The latter 
requires that 7711,2,3 ^ 0.10 eV. For the QD case also one can still ask the question whether 
mi ° r m 3 is the lowest mass, i.e. whether Am\ > or Am\ < 0. 

In the present article we investigate in a model-independent way the possibility of having a 
sizeable value of \U e ^\ as a result of a perturbation of tri-bimaximal neutrino mixing. For 
concreteness, we will use the range of |i7 e 3| in Eq. @ in our analysis. In the light of expected 
deviations from TBM, it represents an interesting and testable benchmark scenario for various 
breaking mechanisms. In this respect, the problem of possible deviations of the neutrino 
mixing matrix from the TBM form has not been studied in detail (see e.g. Ref. [14] for some 
qualitative statements on the subject). Very specific perturbations to TBM in the framework 
of concrete models, with the goal of allowing sizable non-zero #13 ~ 0.1, have recently been 
discussed in Refs. [15]. However, a detailed, quantitative and model-independent analysis, in 
particular in the light of the recent hints for a non-zero U e s, has not been performed before 
and in our opinion is at the present stage both timely and useful. 

More specifically, in this paper we consider values of \U e ^\ ~ 0.1 suggested by Eq. (@J) and 
try to obtain them by starting from TBM. The main challenge is to keep at the same time 
sin 2 #23, and especially sin 2 #12, close to their experimentally determined and thus close to 
the TBM predicted values. As any breaking mechanism introduces correlations between the 
observables, we are able to make characteristic and testable predictions within each case. 
Interestingly, all predictions are connected with CP properties of the lepton sector. 

The paper is organized as follows. In Section [2] we will start by deviating TBM with charged 
lepton corrections and find that CP violation in neutrino oscillations gets constrained to be 
almost maximal by the joint requirement of large \U e %\ and small deviations from sin 2 612 = |. 
Atmospheric mixing deviates from maximal by order sin 2 #23 = \ + C(|^e3| 2 )- Section [3] deals 
with quantum corrections to TBM and shows that only quasi-degenerate neutrinos in the 
Minimal Supersymmetric Standard Model (MSSM) can give rise to sizable \U e s\, while it is 
impossible to produce the required \U e s\ if the effective theory is the standard model (SM). 
Solar neutrino mixing is particularly affected by renormalization effects, but the modification 
of #12 can be suppressed by certain values of the Majorana CP violating phases. These values 
in turn influence the magnitude of the effective Majorana mass in neutrino-less double beta 
((/3/3)oi,-)decay, leading to large cancellations. It is worth noting that the quantum corrections, 
within the context of the MSSM, make sin 2 612 increase, whereas the la range obtained from 
global fits lies below 3. Atmospheric neutrino mixing deviates in general from maximal 
stronger than in the case of charged lepton corrections, namely sin 2 623 = \ + 0(|J7 e 3|). A 
similar but weaker correlation between \U e $\ and the effective Majorana mass in (/3/3)o^-decay 
is found when we explicitly perturb a neutrino mass matrix which without perturbations 
would lead to TBM. This possibility is analyzed in Section [H We also find that in this case 
sizable corrections to sin 2 #23 = \ of order |J7 e 3| are expected. We finally summarize and 
conclude in Section 

2 Breaking Tri-bimaximal Mixing with Charged Lepton Cor- 
rections 

The PMNS matrix is, in general, a product of two unitary matrices, 

U = U\ U v , (6) 
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where U v diagonalizes the neutrino mass matrix and Ui is associated with the diagonalization 
of the charged lepton mass matrix. Several authors have discussed charged lepton corrections 
to various neutrino mixing scenarios [16-20]. It has been shown [18] that, after eliminating 
the unphysical phases, the matrix which diagonalizes the neutrino mass matrix can be written 
as: 

where U u is a "PDG-like" mixing matrix, i.e., 
/ 

U v = 

ii 
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It contains three angles, 8i 2 , #23, and 0\ 3 , and one phase, £. The diagonal matrices P u = 
diag(l, e 1 ^, e tuJ ) and Q u = diag(l, e*°", e"~), in general, cannot be neglected. Note, however, 
that Q v does not affect neutrino oscillation observables [2,21]. The unitary matrix Ug which 
diagonalizes the charged lepton mass matrix can be written as 
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We have used in U,, 
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Let us assume next that U u corresponds to TBM, i.e., U u is given by C/tbm from Eq. ([I]). 
Assume further that the charged lepton corrections are "CKM-like", i.e. that 
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(10) 



with A, B real and of order one. We therefore have in mind here a GUT-like scenario, in 
which tri-bimaximal mixing from the neutrino sector (presumably owing its origin from a 
see-saw mechanism) is corrected by Uj>, which via some quark- lepton symmetry is related to 
the CKM mixing. A natural expectation for A is then that it is kindred to the sine of the 
Cabibbo angle, A ~ sin 0c 1 — 0.227. In scenarios based on SU(5) Grand Unification it often 
happens that a Clebsch-Gordan coefficient of | occurs in between the charged lepton and 
down quark diagonalization, in which case A ~ | sin^c* — 0.076. 

In the case of CKM-like corrections it is straightforward to calculate from U = Uj U v the 



neutrino mixing observables sin 2 #12 
sin #13 = \U e3 



\U e2 \ 2 /(l-\U e3 \ 2 ), sm 2 9 23 = 1^17(1 



u e3 \ 2 



and 



Moreover, it is of interest to obtain the rephasing invariant 
Jcp = Im 



(11) 



which controls the magnitude of CP violation in neutrino oscillations [22], generated by the 
Dirac CP violating phase in the PMNS matrix. In the standard PDG-parametrization of 
the PMNS matrix we have Jcp = g sin 2#i2 sin 2#23 sin 2#i3 cos #13 sin S. The result for the 
observables is 
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plus terms of order A 3 . The magnitude of \U e3 \ is in our analysis fixed by the range in Eq. 
Therefore we can estimate the following interesting range for A: A ~ 0.18lo'o7- 
We can be more general, however, and refine this analysis. To this end, we consider the exact 
and lengthy expression for \U e3 \ and use a random number generator to generate the values 
of \,A,B,u,(ft,ift. We let A vary between and 0.3, the phases between and 2ir, and A, B 
within 0.2 and 5. In order to have a hierarchy in Ut, we take care that sin#f 2 is &t least five 
times as large as sin #23 1 which in turn is at least five times as large as sin 9f 3 . We obtain then 
from the requirement of reproducing the la ranges of the mixing angles given in Eqs. ([31 SJ) 
the range 

0.104 < A < 0.247. (13) 

This is the range for A we will use for the rest of this Section. Interestingly, the sine of 
the Cabibbo angle is included in this range, while one third of it is not. It turns out that 
sin 2 812 can lie anywhere in its currently allowed range given in Eq. ([5]) . In contrast, as 
can be seen in Eq. (|12p . atmospheric mixing receives only small corrections of order A 2 , i.e., 
sin 2 623 = \ + 0(\U e3 \ 2 ). To be quantitative, we find 

0.437 < sin 2 6 23 < 0.533. (14) 

From the expressions for the mixing parameters given in Eq. (I12p . an interesting correlation 
appears [19]: a sizable value of |J7 e 3|> an d therefore of A, introduces a sizable contribution to 
sin 2 #12 of the same order o To be more precise, we have: 

1 2\/2 

- - sin 2 6>i2 ~ — — \U e3 \ cos (p. (15) 

The observed value of | — sin 2 6*12 is at la between 0.002 and 0.039. Thus, cos (ft should lie 
below 0.33, 0.26, or 0.53, if \U e3 \ 2 = 0.016, 0.026, or 0.006. The closer sin 2 12 is to ±, the 
smaller cos (ft is. Consequently, | sin (ft\ is close to one and CP violation in neutrino oscillation 
is "maximal", in the sense that the invariant describing it takes (as a function of |£7 e 3|) almost 
its maximal value El 

We illustrate the phenomenology of this framework in Fig. [TJ The values of the parameters 
in Ut are the same as the ones leading to Eqs. (fl~3l) and flHj). It is easy to see that sin 2 #23 
can have values in a limited interval, and that CP violation is very close to maximal, i.e., 
5 = (ft mod7r ~ 7r/2 or 37r/2. The blue solid lines in Fig. [T] display the maximal value that 
I Jcp| can take. The sign of sin 5 cannot be predicted, because the charged lepton corrections 
to the CP conserving quantities sin 2 9\2 and \U e3 \ fix only cos (ft, whereas CP violation depends 
necessarily on sin (ft. Note that atmospheric neutrino mixing can be maximal. 

Finally, we note an alternative second type of correction from the relation U = Ut U v , namely 
when ut corresponds to TBM and U v is CKM-like [20]. In this case, \U e3 \ ~ sin^/v^- The 
parameter \U e3 \ is therefore governed by the 23-element of U u , which by the same arguments 
as given above for the first case, is expected to be quite small, namely of order A 2 . Even for 
the lowest considered value of \U e3 \ 2 = 0.006 this scenario would require that sin #23 = 0.134, 
a comparably large number, given the GUT-inspired paradigm of "small corrections" in the 
relation u\U v . We will therefore not discuss this possibility further, except for noting two 

2 Similar result holds in the case of hierarchical Ue and U v having bimaximal mixing form [18]. 

3 Solar neutrino mixing then receives correction by the "NLO" term in Eq. (I12|l : sin 2 6*12 = 5 (l + | A 2 ). 
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Figure 1: Charged lepton corrections to tri-bimaximal mixing. The left plot shows \U e ^\ 
against sin 2 #23 (the axes cover the whole la range) while the right plot gives \U e %\ against 
Jcp. The blue solid lines display the maximal value that |Jcp| can take. 



things. First, sizable corrections of order A would arise, in general, for sin #12. Suppressing 
them by choosing a specific value of a CP violating phase is possible, but this phase is not 
related to CP violation in neutrino oscillations. Second, there would also be a very similar 
correlation to Eq. (|15p . namely 

sin 2 023 V2\U e3 \ cos0. (16) 

Note that in this case the atmospheric neutrino mixing angle is correlated with \U e s\ and 
CP violation in neutrino oscillations. We refer to Ref. [20] for more details on this mixing 
scenario. 

3 Breaking Tri-bimaximal Mixing with Quantum Corrections 



Another straightforward breaking mechanism is the application of renormalization group (RG) 
corrections to TBM [19,23-25], which is essential to be considered if the tri-bimaximal scenario 
is assumed to have been generated at some high energy scale. In contrast to charged lepton 
corrections the results now depend on the neutrino mass values and their ordering. In general, 
in the usual PDG-parametrization of the mixing matrix, the corrections to the mixing angles 
can be expressed as [23,24,26]: 

^~^ + a%A r + 0(A 2 ), (17) 

where A is the high scale at which TBM is implemented and A is the low energy scale at 
which measurements take place. We will indicate high scale values by a superscript A in the 
following, and omit for simplicity the superscript A, which would indicate low scale values. 
Hence we have 9± 2 = sin -1 y/T/3, #23 = vr/4 and = 0. We consider the RG evolution of 
the neutrino masses and the mixing parameters in the effective theory and for definiteness 
assume the high scale to be A = 10 12 GeV. The low scale is taken to be A = 10 2 GeV when the 
effective theory is the Standard Model (SM), while we take A = 10 3 GeV when the effective 
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Table 1: Direction of RG correction to the observables sin 2 Q\2 and sin 2 #23 for the SM and 
the MSSM and both possible neutrino mass orderings. 



theory at low energy is the MSSM. The constant C in Eq. (|17p is given by C = —3/2 for 
the SM and C = +1 for the MSSM. The result in Eq. (|17|) is obtained in first order in the 
parameter 
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(18) 



with A e>Al having been neglected since m e>fl <C m T and the vev of the Higgs is taken to be 
u/V2='l74 GeV. 

The dependence on the neutrino mass and mixing parameters is encoded in [23, 24, 27] 



fei2 = — 



y/2 I mi + m 2 e 
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where we have used Q, L j = Ofj, which is correct up to C(A T ), and Am 2 - = m 2 — m 2 . 
running of the masses has been neglected in the above expressions for the kij. The masses 
are decreasing from high to low scale and run as 



m; 



Mi 



(20) 



where Ik is a scalar factor that depends on the SU(2) and U(l) gauge coupling constants 
and the Yukawa matrix in the up quark sector [27-29] and \ii are 0(1) numbers. Thus, 
neglecting the running of masses □ introduces an error 0(A T ) in kij and hence 0(A 2 ) in 6ij. 
One also observes that for |A T | > ((m^ ) 2 — (m 
the C(A T ) terms in the evolution of mj — m\ 



A\2\ 



)/( m o) 1 tli e C(A 2 ) terms dominate over 
[27]. For such cases Eqs. f)19|) will no longer 



be cogent. Thus, for the validity of these equations, we require (mg) 2 A T < (m^) 2 — (m^) 2 , 
which may not be satisfied if (m^) 2 — (m^) 2 is indeed very small. We will therefore use the 
full running equations for the mass matrix itself for the plots and numerical values to be 
presented. Analytical estimates are made with the expressions of the kij and, as we show, 
these estimates can explain the numerical results with a sufficient degree of correctness. 



'Note that the masses appear in both the denominator and numerator of the ki 
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There is a subtle issue involved when we consider k\ 3 in Eq. (|19p . As is seen, fci 3 depends on 
the Dirac CP phase 5 which is unphysical for the case of #13 = at the high scale A. However 
as discussed in [27,30], the value of 5 at this point depends on the values of the masses and 
the Major ana phases and RG evolution takes care of that automatically. 
For analytical estimates, it is convenient to consider the shift of the mixing angles 0y from 
their initial values. From the above expressions for the kij, and in the limit of \kij A T | <C 1, 
one obtains the following expressions for the observables: 

|sm0 13 |~|C&i 3 A T | , sm 2 9 23 ~^-Ck 23 A T , sin 2 12 - i ~ ^ C k 12 A T . (21) 

In the spirit of our analysis we require (see Eqs. ([3l H|)) that \C fc 13 A T \ = 0.077 — 0.161, while 
-C k 12 A T = 2.8 • 10~ 3 - 4.2 • 10~ 2 . Note that for the la range we are taking, C k 12 A T (and 
therefore C) is supposed to be negative. Hence, within the MSSM the required deviation 
from TBM cannot be realized. Therefore we use the 3a ranges for sin 2 6± 2 and sin 2 #23 f° r 
the purposes of illustration. If indeed the trend of sin 2 6*12 < | continues then it will not 
be possible to account for a high scale value of sin 2 9\ 2 = | solely by RG effects within the 
MSSM0. 

From Eqs. (|17p and (119p it is evident (and well-known) that whether the angles Oij will decrease 
or increase during evolution will depend on the effective theory (SM or MSSM, through the 
factor C) and also on the sign of Am 31 (~ Am 2 2 ) for #23. Table Q] summarizes the direction 
of the correction for the observables sin 2 8i 2 and sin 2 #23- Since k\ 2 is always positive, #12 at 
the low scale is always larger (smaller) than that at high scale for the MSSM (SM). The size 
of the RG corrections will depend on the values of the Majorana phases, the neutrino masses 
and, in case of the MSSM, tan /3. 

One can at the outset make some interesting observations from Eqs. (I19j) . It is easy to 
see that only quasi-degenerate neutrinos will be able to lead to values of \U e s\ around 0.1. 
Note also that in this case the running of solar neutrino mixing is in general enhanced by 
a factor \Am\\/ Am^ with respect to the running of the other mixing angles. We will see 
that suppressing the running of sin 2 0\ 2 and hence of the factor k\ 2 by suitable values of the 
Majorana phase a 2 around ir has interesting phenomenological consequences in (/3/3)o;,-decay. 
We further can expect that the deviation from maximal #23 is of the same order than the 
deviation from zero \U e 3\. In the following we will quantify these statements. 
We have performed a detailed analysis, by numerically solving the RG running equations for 
the effective neutrino mass matrix and then diagonalizing it to extract the masses, mixing 
angles and phases at low scale. At the high scale A, the angles 9fj are fixed by the requirement 
of the TBM scenario, while the masses and the CP phases are chosen randomly so that after 
the RG evolution at low scale the parameters are consistent with the chosen ranges of the 
current experimental data. We have used the following ranges for the high scale values of the 
mass-squared differences: (Arag) A = 10 -6 — 10 -3 eV 2 and |Am^| A = 1.5 x 10~ 3 — 10 -2 eV 2 , 
while the phases are varied over the full range of — 2tt. 

Starting with the SM, Fig. [2] shows the allowed region in the tuq - sin 2 #13 plane at the low 
scale A, after performing the RG evolution, for both the normal (left panel) and inverted 
(right panel) mass orderings. Recall that m 2 , ^ lAr/r^J is the common neutrino mass scale 

5 Note that for the same reason any initial value of sin 2 612 > ^ (including bimaximal mixing) at high scale 
is excluded unless of course highly model-dependent see-saw threshold effects [30,31] are taken into account. 
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Figure 2: The running of | t/e3 1 2 = sin 2 #13 in SM for both the normal (left panel) and the 
inverted (right panel) mass orderings. The high scale values of mixing angles are kept fixed at 
TBM values while the masses and phases are varied randomly such that after RG evolution 
the parameter values are within current experimental ranges. 



for quasi-degenerate neutrinos. As can be seen, to generate values of \U e s\ within the range 
of interest, neutrino masses should exceed the direct limit of 2.3 eV from tritium decay [32], 
and hence also the more stringent but model-dependent limits from cosmology. We conclude 
that a high scale value of #13 = is incompatible with the indicated range of |J7 e 3|. The 
dependence of this statement on the initial values of and #23 is moderate and hence this 
statement is valid in general. 

We will focus on the MSSM in what follows. As already stated above, we require the 3a 
ranges of the oscillation parameters to be satisfied, because, strictly speaking, the MSSM 
cannot reproduce the la range, due to its prediction of sin 2 #12 > 3-. 

Fig. [3] shows the allowed region in the uiq - sin 2 #13 plane, when the effective theory is the 
MSSM, for tan [3 = 5, 20 and the normal mass ordering. The left panel shows that sin 2 #13 
lies in the required range when 0.8 eV ^ mo ^ 1.4 eV for tan/3 = 5, while the allowed 
mass range becomes 0.2 eV ^ niQ ^ 0.34 eV for tan/3 = 20, as can be seen from the right 
panel. Thus, the relevant range of mo tan (3 is given by (see below for analytical estimates) 
4.1 ^ (mo/eV) tan/3 ^ 6.9. Hence the allowed mass ranges depend strongly on tan/3 and for 
higher values of tan /?, lower values of mo are sufficient to produce the required running of 
#13. It has been checked that for a fixed tan (3 value, there is no significant dependence on the 
mass ordering, other than the direction of the correction to #23- From the allowed mass ranges 
obtained in Fig. [3] it is seen that to have sin 2 #13 in the la range under consideration, we need 
the neutrinos to be quasi-degenerate even for the MSSM with tan /3 = 20. Fig. |4]shows scatter 
plots of the allowed region of the neutrino mass scale mo and the Major ana phase 02, which 
is particularly important for the running of #12 [27] (see also [33]). We compare the allowed 
regions at high and low scale for a normal mass ordering and tan (3 = 5, 20. The scattered 
plots obtained for the inverted mass ordering show the same characteristics. We see that [02! 
is restricted in a narrow region around [c^l = tt for all cases. 

In order to explain the plots analytically we consider Eqs. (119p in the QD regime tuq 3> Am\ 
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Figure 3: Scatter plots showing the running of sin 2 #13 with mo, for MSSM with normal mass 
ordering and tan (3 = 5, 20. The high scale mixing angles are fixed at TBM and the masses 
and phases are varied randomly such that after RG evolution the parameter values are within 
current experimental ranges. For a given tan (3, the allowed regions are the same as above for 
the inverted mass ordering. 



to obtain: 



,2 



sin 2 0i 2 --) ~ -CA T (1 + 00802)^3^-, (22) 
\/2 m 2 

I sin6>i 3 |QD — — C A T ° 9 |(1 + it!) cos(£ + a 3 - a 2 ) - cos(5 + a 3 ) + cos<5| , (23) 
3 Am^ 

where R = Attiq/ Am\. From Eq. (|22p one can understand that the low energy constraint 
on sin 2 9\2 from the current experimental data restricts | ct2 1 to remain close to it, as shown in 
Fig. [J making (1 + cos a%) small so that there is less running of #12 even with large neutrino 
masses. The plots in Fig. [J] further show that is also close to it and that «2 stays close to 
it in the course of its RG evolution. This can be estimated from the fact that the running of 
«2 can be expressed as a\ ~ + 02 A T [24, 27] with 02 — —2/(3 AttIq) m^ m^ sin a%. From 
Eq. (|22p we note that the maximum running for 612 is obtained for 02 = 0. In absence of any 
lower bound on #13 this value was still allowed [24]. However if we put «2 = in Eq. (|23p then 
the running of #13 is suppressed by the factor \R\ = Am^ ) /\Am\\. Thus, the requirement of 
large running of #13 disfavors 02 = and further strengthens the bound in the 0.2 - uiq plane. 

In the limit of 02 = it and quasi-degenerate neutrinos, the maximum value of | sin 9\^\ that 
can be achieved starting from 6^ = can be estimated from Eq. (|23p as 



2V2 



-2 



I sin 13 | QD < ±Xf C A T -pL- (1 + R) > (24) 

Am. 



with 03 = 0, 5 = ±7r or 03 = ±7r, 5 = 0. Thus, from Eq. (|24|) one can estimate that 
013 > 0.077 requires m £ 2.66 eV for the SM and m Z 0.72 (0.18) eV for tan/3 = 5(20) 
with the MSSM. The estimates are in good agreement with the allowed mass ranges obtained 
from Fig. [2] and Fig. El respectively. 
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Figure 4: Scatter plots in the mo - \ot2\ plane for MSSM (tan/3 = 5,20) and normal mass 
ordering, both for high (black circles) and low (red squares) energy scales. The high scale 
mixing angles are fixed at the TBM values and the masses and phases are varied randomly at 
the high scale so that the low energy parameters are consistent with the current experimental 
data. The data for the inverted mass ordering shows the same variation. 



Fig. [5] shows the correlation between the low scale values of sin 2 #13 and sin 2 #23. For normal 
ordering sin 2 623 > \ , whereas for inverted ordering sin 2 623 < \ ■ For normal ordering #13 and 
#23 are correlated, i.e., a higher value of #13 requires a higher value of 623- For the inverted 
ordering the predicted values of the two angles are anti-correlated. The plots obtained with 
tan (3 = 20 are identical to those shown in Fig. [5] for tan (3 = 5, when the mass ordering is the 
same. For a different tan (3 the value of mo adjusts itself to comply with the low energy cuts 
on the parameters and the allowed points in the sin 2 #23 — sin 2 #13 plane remain same. 
We note here that maximal atmospheric neutrino mixing is not possible. To be more quanti- 
tative, we find that 



(25) 



0.55 < sin 2 #23 < 0.64 for Am| x > , 
0.33 < sin 2 #23 < 0.45 for Am^ < , 

independent on the value of tan (3. 

In Fig. [6] we plot the effective Majorana mass 

(m) = cos 2 #13 mi cos 2 #12 + m 2 sin 2 9 12 e ia2 + m 3 tan 2 # 13 e i{a3+2S) , (26) 

which governs the rate of (/3/3)o^-decay at low energy. The scatter points show the values of 
(m) allowed by the low energy neutrino oscillation data after RG analysis. The solid (black) 
lines indicate the maximum and minimum possible values of (m) at low scale for a given 
mo, obtained by varying the oscillation parameters in their current 3<r range and the phases 
between to 27r. The plots show that the effective mass obtained after RG analysis lies close 
to its minimum allowed range. As can also be seen from Fig. [6l for tan (3 = 5, (m) takes 
values between 0.26 and 0.50 eV, to be compared with the general upper and lower limits 
of 0.2 eV and 1.4 eV. If tan/? = 20, then 0.07 eV < (m) < 0.11 eV, while in general the 
effective Majorana mass could be in between 0.05 eV and 0.34 eV. As can be estimated from 
Eq. (|26p . the maximum value that (m) can achieve for quasi-degenerate neutrinos is when 



11 





0.026 r 




0.022 








0.018 








.B 




GQ 


0.014 




0.01 




0.006 




sis 



0.026 



0.022 



0.018 



0.006 



Inverted 
tan/3 = 5,20 



m o.oi4 v. j*; 



0.01 



0.33 0.35 0.37 0.39 0.41 0.43 0.45 
,2, 



sin Q< 



Figure 5: Scatter plots showing the correlation between sin 2 #13 and sin 2 #23 for MSSM with 
tan/3 = 5, 20. The left panel shows the case with normal mass ordering, while the right panel 
is for the inverted mass ordering. The high scale mixing angles are fixed at TBM and the 
masses and phases are varied randomly such that after RG evolution the parameter values 
are within current experimental ranges. 



ot2 = «3 + 25 = 0, and is given by [34] mo, whereas the minimal value is obtained when 
ct2 = «3 + 25 = 7r: 



(m)™ n ~ m (cos 20 12 - 2 \U e3 \ 2 /(l + tan 2 9 12 ] 



(27) 



As we have seen, RG evolution combined with low energy constraints imply QD neutrinos 
with «2 close to tt. In this limit 



( m )ob n - m o cos2 #13 



cos 26»i2 + tan 2 6 13 e i( - as+25) 



(28) 



and since #13 is small at all energy scales, expanding in powers of |J7 e 3| = sin #13 we can 
write [34] 



MqiT-^o (cos20 12 + C>(|Z7 ( 



e3 



(29) 



Thus, RG evolution constrains (m) towards the minimum allowed value, which is confirmed 
by the figure. 



We can give very simple forms of the neutrino mass matrix in the flavor basis satisfying the 
above constraints. In general the mass matrix generating TBM reads 



I A B B \ 

■ \{A + B + D) \{A + B-D) . (30) 

V • • \{A + B + D) ) 

The parameters A,B,D are in general complex and functions of the neutrino masses and 
Majorana phases: 



A = \ (2mi+m 2 e~ ia2 



B = g \ m-2 ( 



-totl 



mi ) , D = m?, e 



-ia 3 



(31) 
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Figure 6: Scatter plot for the effective neutrino mass (to) that contributes to neutrino-less 
double beta decay as a function of too, in MSSM with tan (3 = 5, 20 and normal mass ordering. 
The solid (black) lines indicate the maximum and minimum possible values of (m) for given 
mo, obtained by varying the oscillation parameters in their current 3a range and the phases 
between to 2tt. The cases with inverted ordering show identical characteristics. 



Note that the sum of the elements in each row, and in each column, equals A + 2B = 771,2 e~ ia2 . 
Now to estimate the texture of to„ at high scale let us insert 777^2,3 = m o, TBM and a 2 = ^- 
It follows 



_3_ 

777 



1 -2 

• i(-l + 3e" ia3 ) 



-i(l + 3e" iQ3 ; 
\ (-1 + 3e~ ia3 ) 



1 



I V 



-2 -2 
1 -2 
1 

-2 -2 
-2 1 
• -2 



for «3 = . 



for q: 3 = 7T , 



where we have set two specific values of a^. Corrections to these expressions are of order 



At77q/777o, J Am a /too and hence small for QD neutrinos. 



4 Breaking Tri-bimaximal Mixing explicitly 

We can explicitly break TBM by perturbing the neutrino mass matrix. As for RG effects, we 
will see that there is crucial dependence on the neutrino mass ordering and values of neutrino 
masses. In its general form, m u leading to TBM is given in Eq. (|30p . A possible strategy to 
perturb TBM, outlined in detail in Ref. [35], is to modify the mass matrix in the following 
way: 

/ A(l + £l ) B(l + e 2 ) 



\{A + B + D)(l + e i ) 



B(l + e 3 ) \ 
i( A + B-D)(l + e 5 ) . (32) 
V • • ±(A + B + D)(l + e 6 ) J 

The complex perturbation parameters ej are taken to be |ej| < 0.2 for i = 1 — 6 with their 
phases 4>i allowed to lie between zero and 2tt. In case of a normal hierarchy, one finds [35] that 
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Figure 7: Scatter plot for an explicitly broken TBM mass matrix of | U e 3 | against the smallest 
neutrino mass for the normal (left) and inverted (right) mass ordering. 



[ t/e3 1 2 is of order e 2 R, where e is the magnitude of one of the e^, and R = Amg/Am^. Hence, 
a too small value of |C/ e 3| 2 is generated in this case. It turns out that at least mi ~ 0.015 eV 
is required in order to generate |C/ e 3| 2 above 0.006. This is illustrated in Fig. 
Such values correspond to a scenario with a partial mass hierarchy: mi ~ m.2 ^ rn-^. With the 
increase of m\ starting from 0.015 eV, the maximal value of \U e ^\ grows almost linearly with 
m\. In contrast, in the case of inverted hierarchy (ordering), one can generate large values 
of \U e z\ even for a vanishing value of the smallest neutrinos mass 7713. For quasi-degenerate 
neutrinos, obviously, sizeable values of \U e 3\ ~ 0.1 can also be generated. In addition, in 
the cases of neutrino mass spectrum with partial hierarchy, with inverted hierarchy and of 
quasi-degenerate type, there exists a correlation between the effective Majorana mass in 
(/3/3)oi/-decay and the value of |t7 e 3| thus generated. 

To illustrate the above comments, consider the following analytic estimates in the case of 
spectrum with inverted ordering. We first set 02 = vr. In this case one has A ~ 



and B ~ — 2wAmj^/3. Consider now a perturbation of the form 



[A B(l + e) B(l-e) \ 

■ \{A + B + D) \{A + B-D) 



\ 



(33) 



UA + B + D) J 



for real e, either negative or positive. In this case we get the largest effects on |f7 e 3| and 
sin 2 23 [35]: 



\U, 



e3 



8 16 m.3 
81 + 27 



< 10" 2 and 



sin 



723 



~^e 2 >0.18. (34) 
9 



Note that A 



~ 0.016 eV is the minimal possible value of the ( / 9/3)o ! y-decay 



effective Majorana mass in the case of spectrum with inverted hierarchy under discussion. 
The cancellation arises due to the chosen CP conserving value of the Majorana phase «2- 
In the other extreme case of 02 

at most of order (e.B/^4) 2 ~ 10 _b and therefore completely negligible. 



|2 : r 



0, we have B/A ~ g Attiq/ Am\ and we find that |£7 e 3r |S 
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Figure 8: Scatter plot of \U e ^\ against sin 2 #23 as well as of |?7 e 3| against (m) for an explicitly 
broken TBM mass matrix in case of an inverted hierarchy. The left plots show the cases 
«2 = and «2 = 7T) the right plots have free ai- Indicated are also the lcr ranges of the 
oscillation parameters, and the upper and lower limits of the effective mass. 



In the case of perturbed and rr entries of m u , 



m. 



( 



\ 



A 



B 



UA + B + D){l-e) 



B 

\{A + B-D) 
\{A + B + D){l + e) j 



(35) 



the largest possible deviation of #23 from tt/4 is obtained for «2 = 0: | sin 2 #23 — 7j| ^ e/2 ~ 0.1. 
We conclude that [35], if initially the phase «2 takes a CP conserving value of n and for an 
inverted hierarchy neutrino mass spectrum, perturbed TBM leads to values of the (/3/3)o- 

decay effective Majorana mass close to the minimal one, (m) = c 2 3 J Am\ cos 29\2- These 
values are correlated with sizable values of \U e s\ and relatively large deviations from maximal 
atmospheric neutrino mixing. The benchmark value of \U e 3\ from Eq. can be reconciled 
with minimal allowed values of the effective Majorana mass. In contrast, if the effective 

Majorana mass (m) is close to its possible maximal value, {m) ~ cf 3 \J Am^, negligible values 
of |£/ e 3| are predicted. Hence, the benchmark value of \U e 3\ from Eq. (|4|) cannot be reconciled 
with values of (m) close to its maximal value. The expected deviation from sin 2 #23 = \ is 
also smaller than in the previous case. It turns out, however, that the case of free «2 7^ or 
7r allows non-minimal values of (m) for sizeable \U e ^\ as well (see below), which means that 
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the correlations discussed above rely on extreme initial values of a.2- 

In Fig. [8] we show scatter plots resulting from a corresponding numerical analysis. We have 
diagonalized Eq. (|32p , where we have taken random values for the complex a , by starting with 
m2 = 0.051 eV, mi = 0.0502424 eV and 777.3 = 0.01 eV. We required the resulting oscillation 
observables to lie in their 3<r ranges. We have also chosen as initial values 02 = 0, ir, but let 
02 vary freely as well. The largest and smallest possible values of the effective mass (m) are 
approximately 0.059 eV and 0.0135 eV, respectively, and we have indicated them in the figure. 
The analytical estimates from above are confirmed here. Fig. [9] shows the same analysis for 
quasi-degenerate neutrinos, where we have started with 7773 = 0.10 eV, 7772 = 0.08778 eV and 
mi = 0.08735 eV. The effective Majorana mass (m) in this case lies between 0.023 eV and 
0.105 eV. In both cases it is evident that the value of #23 is not a good discriminator. In 
particular, maximal atmospheric neutrino mixing is always possible. 



5 Conclusions and Summary 

Tri-bimaximal mixing provides a very close description of neutrino mixing angles. However, 
the present hint of non-zero #13 coming from analyses of the global neutrino oscillation data 
may indicate that it is broken. In this paper we consider three breaking mechanisms from 
exact TBM - charged lepton corrections, radiative corrections and explicit breaking. While 
the deviation from maximal sin 2 623 = \ is allowed by the data to be of the same order 0.1 
as the values of \U^\ that we study, the challenge is to simultaneously keep the deviations 



16 





charged leptons 


renormalization (MSSM) 


explicit breaking 


sin 2 (9 2 3 


0.44 - 0.53 


0.55 - 0.64 (Ami > °) 
0.33 - 0.45 (Am| < °) 




Peal 


\ 

" V2 


oc -^\(l + tan 2 /3) 
Ami 


oc e (IH) 
ocemi/^Ami (PD/QD) 


mass 




QD: m tan^ ~ (4 - 7) eV 


IH, PD, QD 


(m) 




mo c 2 3 cos 29\2 


mo c 2 3 cos 2#i2 (QD) 
^/Ami cf 3 cos 26»i2 (IH) 


CP 


oscillations: almost 
maximal CP violation 


«2 — 7T 


large \U e ^\ requires 
suppressed (m) only 
when initially cti ~ 7r 



Table 2: Requirements on and predictions of the three breaking scenarios in order to generate 
the la range 0.077 < |Z7 e 3| < 0.161. IH denotes inverted hierarchy, while PD stands for a 
partial hierarchal and QD for a quasi-degenerate mass scheme. 



from sin 2 6\2 = \ of order |?7 e 3| 2 or below. For definiteness we choose in our analysis the la 
range given in Eq. 0.077 < |J7 e 3| < 0.161. The main results of this work are summarized 
in Tabled 

We first assume CKM-like charged lepton corrections from Ug to U u corresponding to tri- 
bimaximal mixing. The correction parameter A, which is the sine of the 12-rotation in the 
usual parametrization of Ug, can be restricted as 0.104 < A < 0.247 from the current la ranges 
of the mixing angles. We note that the sine of the Cabibbo angle is included in this range, but 
one third of it is not. In this picture pes] — A/v2- A sizable value of U e ^, therefore implies 
a sizable A. Suppressing the leading (0(A)) correction to sin 2 #12 is possible by choosing the 
Dirac CP phase in neutrino oscillations to be ir/2 or 37r/2 corresponding to maximal CP 
violation in neutrino oscillations. The charged lepton corrections to tri-bimaximal mixing do 
not depend on the neutrino mass values and their ordering. The atmospheric neutrino mixing 
parameter sin 2 #23 is deviated from | by terms of order |C/ e 3| 2 . To be precise, it is within the 
range 0.44 ^ sin 2 #23 ^ 0.53. In particular it is allowed to be maximal. 

Generating a large \U e s\ via radiative corrections implies quasi-degenerate neutrinos, namely 
mo 2.6 eV for the SM when U e ^ = at the high scale. Thus, the current neutrino mass 
limits rule out a possible RG origin of p&\ — 0.1 in the SM. For the MSSM one requires mo ^ 
0.8 (0.2) eV with tan/3 = 5 (20). The implied constraints of the la range of |£7 e 3| on mo and 
tan (5 can be summarized as 4 ^ (mo/eV) tan (3^7. The running in the MSSM predicts that 
sin 2 6*12 increases from its initial high scale value. Large running of 6*13 to generate \U e 2,\ — 0.1 
together with the requirement that sin 2 #12 is within its current 3a range forces the Majorana 
phase «2 — 7T- Interesting correlations are also obtained between sin 2 #13 and sin 2 6*23. The 
latter parameter is necessarily non-maximal and lies in the range 0.55 < sin 2 #23 < 0.64 for a 
normal ordering and 0.33 < sin 2 #23 < 0.45 for an inverted ordering. The RG evolved effective 
neutrino Majorana mass observed in neutrino-less double beta decay is found to lie close to 
its minimum allowed value because of the constraint of ai — tt. In case of tan/? = 5, (m) lies 
between 0.26 and 0.50 eV, to be compared with its general upper and lower limits of 0.2 eV 
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and 1.4 eV. If tan/3 = 20, then 0.07 eV (m) ^ 0.11 eV, while in general the effective mass 
could be in between 0.05 eV and 0.34 eV. 

We also consider the possibility of deviating from tri-bimaximal mixing by adding explicit 
breaking terms to the neutrino mass matrix, i.e., every entry is multiplied with an individ- 
ual factor 1 + 6j. For this mechanism to generate sizable \U e s\ the neutrino mass spectrum 
has to be partially degenerate, or quasi-degenerate, or with inverted hierarchy. Atmospheric 
neutrino mixing is allowed to take any of its currently allowed values, including sin 2 #23 = \. 
In this breaking scenario the requisite sizeable \U e ^\ value cannot be reconciled with initial 
maximal values of the effective Majorana mass governing neutrino- less double beta decay, 
corresponding to the three indicated types of neutrino mass spectrum. 

To sum up, the CP violating phases in the neutrino mixing matrix play a crucial role for having 
only relatively small corrections to #12 when large corrections to U e ^ = are generated. This 
interesting fact together with the predictions for #23 may be used to distinguish breaking 
scenarios to tri-bimaximal mixing. 
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